Supersolidity in a Bose-Holstein model 



o 
o 

(N 



Sanjoy Datta and Sudhakar Yarlagadda 
CAMCS and TCMP Div., Saha Institute of Nuclear Physics, 1/AF Salt Lake, Kolkata-64, India 

(Dated: July 21, 2009) 

We derive an effective d-dimensional Hamiltonian for a system of hard-core-bosons coupled to 
optical phonons in a lattice. At non-half-fillings, a superfluid-supersolid transition occurs at inter- 
mediate boson-phonon couplings, while at strong-couplings the system phase separates. We demon- 
strate explicitly that the presence of next-nearest-neighbor hopping and nearest-neighbor repulsion 
leads to supersolidity. Thus we present a microscopic mechanism for the homogeneous coexistence 
of charge-density-wave and superfiuid orders. 
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Introduction: Coexistence of charge-density-wave 
(CDW) and superconductivity is manifested in a variety 
of systems such as the perovskite type bismuthates (i.e., 
BaBi03 doped with K or Pb) [1], quasi- two-dimensional 
layered dichalcogenides (e.g., 2H — TaSe 2 , 2H — TaS2, 
and 2H-NbSe 2 ) 0, etc. What is special about the 
above systems is that they defy the usual expectation 
that the competing CDW and superconductivity orders 
occur mutually exclusively. CDW and superconductiv- 
ity are examples of diagonal long range order (DLRO) 
and off-diagonal long range order (ODLRO) respectively. 
Here, DLRO breaks a continuous translational invariance 
into a discrete translational symmetry, whereas ODLRO 
breaks a global U(l) phase rotational invariance 

Another interesting example of DLRO-ODLRO con- 
currence is a supersolid (SS). A SS state is characterized 
by the homogeneous coexistence of two seemingly mu- 
tually contradictory phases, namely, a crystalline solid 
and a superfiuid (SF). Here all particles simultaneously 
participate in both types of long range order. It has 
been conjectured long ago that solid helium-4 may ex- 
hibit supersolidity due to Bose condensation of vacancies 
present 0, Only recently, Chan and Kim (using a 
torsional oscillator) observed a decoupling of a small per- 
centage of the helium solid from a container's walls Q. 
They interpreted this as supersolidity. This discovery led 
to studies of bosonic models in different kinds of lattice 
structures [3] and with various types of particle interac- 
tions Phenomenologically, Ginzburg-Landau theory 
[§] and a quantum solid based on Gross-Pitaevskii equa- 
tion [H| have been used to study supersolidity. 

In this paper, we address the above ongoing puzzles 
by studying the quantum phase transitions exhibited by 
hard-core-bosons (HCB) coupled to optical phonons. To 
this end, in contrast to the above treatments, we em- 
ploy a microscopic approach involving a minimal Bose- 
Holstein (BH) lattice model. Examples of real systems 
describable by our BH model are as follows. In the bis- 
muthate systems, the observed valence skipping of the 
bismuth ion is explained by invoking non-linear screen- 
ing which is said to produce a large attractive interaction 
resulting in the formation of local pairs or HCB 11 , . 



Such HCB couple to the cooperative breathing mode of 
the oxygen octahedra surrounding the Bismuth ion. In 
a helium-4 crystal, vacancies produce a local distortion 
and can be treated as HCB coupled to finite frequency 
phonons. Furthermore, the concentration of vacancies is 
very small Q and hence direct interactions among these 
HCB is negligible. Lastly, for dichalcogenides such as 
NbSe2, where homogeneous coexistence of the two long 
range orders has been unambiguously established [l3j ]. 
our results on the BH model should be quite relevant. 

Starting with our BH model, we derive an effective d- 
dimensional Hamiltonian for HCB by using a transparent 
non-perturbative technique. The region of validity of our 
effective Hamiltonian is governed by the small parame- 
ter ratio of the adiabaticity t/u>o and the boson-phonon 
(b-p) coupling g. The most interesting feature of this 
effective Hamiltonian is that it contains an additional 
next-nearest-neighbor (NNN) hopping compared to the 
Heisenberg xxz-model involving only nearest-neighbor 
(NN) hopping and NN repulsion 3, HI ■ We employed a 
modified Lanczos algorithm [l|| (on lattices of sizes 4x4, 
y/18 x a/18, V20 x \/20, and 4x6) and found that (ex- 
cept for the extreme anti-adiabatic limit) the BH model 
shows supersolidity at intermediate b-p coupling strengths 
whereas the xxz-model produces only a phase separated 
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FIG. 1: Quantum phase diagram at various particle numbers 
N p , t/ujo = 1-0 and (a) NNN hopping J2 > in our effective 
BH Hamiltonian [Eq. ©] and (b) J2 = (xxz-model). Su- 
persolidity occurs only in (a). The continuous lines in (a) and 
(b) are guides for the eye. 
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(PS) state. Here we present the calculations for only a 
4x4 lattice (see Fig. [T]) . Our main results are depicted in 
Fig. [1] (a), where supersolidity manifested by local pairs 
(as HCB) implies homogeneous coexistence of CDW and 
superconductivity. 

Effective Bose-Holstein Hamiltonian: We start 
with a system of spinless HCB coupled with optical 
phonons on a square lattice. This system is described 
by a BH Hamiltonian [l7| 

H = -t y] b]bj + s+wo }] of fflj+gctfr rij (oj+fflj), (1) 

3,8 j j 

where 6 corresponds to nearest- neighbors, ujq is the op- 
tical phonon frequency, bj (a,-) is the destruction op- 
erator for HCB (phonons), and rij = bjbj. Then we 
perform the Lang-Firsov (LF) transformation [HI 
on this Hamiltonian and this produces displaced simple 
harmonic oscillators and dresses the hopping particles 
with phonons. Under the LF transformation given by 



bj and 



e s He- s = H + H' with S = -gJ2i n ii a i ~ 4) 
dj transform like fermions and phonons in the Holstein 
model. This is due to the unique (anti-) commutation 
properties of HCB given by 

[Mi] = M]]=0, iorijtj, 

{h,b\} = 1- (2) 
Next, we take the unperturbed Hamiltonian Hq to be [HI 

H = -J\ 2J b]b j+ s + E a ) a j ~ 9 2uj o E n i> ( 3 ) 

3,8 3 3 

and the perturbation H' to be 

H' = -J^b]b j+s {S^Si-l}, (4) 

3\8 

where S± — exp[±g(aj — a,j+s)], Ji — te~xjp(—g 2 ), and 
g 2 u}Q is the polaronic binding energy. We then follow the 
same steps as in Ref. [lU to get the following effective 
Hamiltonian in d-dimensions for our BH model 

H e = -g 2 Ul ^ n j - Jl E b ] b 3+ s 
3 3,8 

3,8,8'^8 



u j+5 ,b j+ s-0.5J z 



£%(1 -n j+s ), (5) 
3,8 



where J z = ( Jf /u )[Af x (g) + 2f 2 (g)} and J 2 = 
{JlMHg) with h(g) = EZi9 2n /(n\n) and f 2 (g) = 
12ri=i Em=i g 2 ( n+m ) /[n\m\{n + m)j. Here we would like 
to point out that, as shown in Ref. 12CJ, the small param- 
eter for our perturbation theory is t/(gu3o). 

Long range orders: Diagonal long range order 
(DLRO) can be characterized by the structure factor de- 
fined in terms of the particle density operators as follows: 



Off-diagonal long range order (ODLRO) in a Bose- 
Einstein condensate (BEC), as introduced by Penrose 
and Onsager Q], is characterized by the order parame- 
ter (bo) = y/ (no}e l6 where n is the occupation number 
for the k = momentum state or the BEC. ft is useful 
to define the general one-particle density matrix 



P(i,j) 



^■>-^E e(k - R ^ q - Rj) 

k,q 



(7) 



where () denotes ensemble average. Eq. ([7]) gives the 
BEC fraction as 



n b 



[ n o) 



E 

1,3 



(8) 



fn general, to find rij,, one constructs the generalized one- 
particle density matrix p and then diagonalizes it to find 
out the largest eigenvalue. To characterize a SF, an im- 
portant quantity is the SF fraction n s which is calculated 
as follows. Spatial variation in the phase of the SF order 
parameter will increase the free energy of the system. We 
use a linear phase variation 0(x) = 9q(x/L) with 9q being 
a small angle and L the linear dimension in x-direction. 
This is done by imposing twisted boundary conditions 
(TBC) on the many-particle wave function. At T = K, 
we can write the change in energy to be 



E[0 O ] - E[0] = ImNj 



m 



Then the SF fraction is given by [2JJ, |2 
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N„t 



P l eff 



E[9 ] - E[0] 



(9) 



(10) 



where t e ff = h 2 /2m. For our Hamiltonian in Eq. (fSJ), we 
find t e ff = Ji + 8J2. The phase variation (taken to be 
the same for BEC and SF order parameters) is introduced 
in our calculations by modifying the hopping terms with 
bj — ► bj ; exp[ix • Rj9 Q /L] which is gauge-equivalent to 
TBC. 

Results and discussion: We employ the mean field 
analysis (MFA) of Robaszkiewicz et al. [lj] to study the 
phase transitions dictated by the effective Hamiltonian of 
Eq. |5]). We obtain the following expression for the SF- 
PS (SF-CDW) phase boundary at non-half-filling (half- 
filling): 



Jz_ _ 3J2 _ 1 + (2n - 1) 
2 Ji Ji ~ 1 - (2n - 1) 



2 ■ 



(11) 



Eq. (jTTJ) leads to the same phase diagram as that (for the 
xxz-model) in Ref. |l4j but with J^^ = J z /2 — 3J 2 as 
the y-coordinate instead of J^f = J z /2. Thus, within 
mean field, NNN hopping does not change the qualitative 
features of the phase diagram; it only increases the criti- 
cal value of J z j J l at which the transition from SF state 
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FIG. 2: Free energy at different fillings and different values 
of g for (a) Ji > 0, J 2 / 0, and t/cu = 0.1; and (b) Ji = 0, 
J 2 / 0, and t/u> = 1-0. 



to PS or CDW state occurs. However, as demonstrated 
by the numerics below at non-half-fillings, MFA fails to 
capture the supersolid phase. 

We studied the stability of the phases by examining the 
nature of the free energy versus N p curves at different b-p 
couplings g (see Figs. [2] and El and by using an analysis 
equivalent to the Maxwell construction. For the system 
at a given N p , g, and t/ui^, if the free energy point P 
lies above the straight line joining the nearest two stable 
points Q and R (lying on either side of P) on the same 
free energy curve, then the system at P breaks up into 
two phases corresponding to points Q and R. 

For the range of parameters that we considered (i.e., 
0.1 < t/ujQ < 1 and g > 1), the small parameter 
t/(guJo) < 1. The behavior of the system for i/w = 0.1 is 
the same for both J 2 = and ^ Jij J\ [= (Ji/^o)fi{g)] 
because J2/J1 is negligible in the latter case. Further- 
more, for J2 = 0, the behavior of the system as a func- 
tion of g is qualitatively the same for all values of t/cuo as 
there is only one dimensionless parameter J z / Ji involved 
in Eq. ©. 

We will now analyze together, in one plot, the 
quantities nb, n s , and the normalized structure factor 
S*(tt,it) = S(n,TT)/S max {Tr,TT) where S max (tt, tt) corre- 
sponds to all particles occupying only one sub-lattice. 





FIG. 3: Plot of free energy for different number of particles 
and various values of g when Ji 7^ 0, J 2 7^ 0, and t/u>o = 1.0. 



FIG. 4: Comparative plots of S*(ir, n), rib, and n s when Ji 7^ 
0, J 2 =£ 0, and (a) t/u = 0.1, N p = 8; (b) t/w Q = 1.0, N p = 8; 
(c) t/cuo = 0.1, N p = 5; and (d) t/uj = 1.0, N p = 5. 



At half-filling, we first observe that the system is either 
a pure CDW or a pure SF. For a half-filled system (i.e., 
N p = 8) at J 2 ^ and t/u = 0.1 (t/cj = 1.0), we 
can see from Fig. 2] (a) [Fig. 0] (b)] that the system un- 
dergoes a sharp (first order) transition to an insulating 
CDW state at g c sa 2.15 (g c s» 1.55). At g = g c , while 
there is a sharp rise in iS*(7r,7r), there is also a concomi- 
tant sharp drop in both the condensation fraction rib and 
the SF fraction n s . Furthermore, while n s actually goes 
to zero, rib remains finite [as follows from Eq. ©] at a 
value 1/N = 1/16 which is an artifact of the finiteness 
of the system. Larger values of t/ojo for a half-filled sys- 
tem leads to lower values of g c . This is in accordance 
with the MFA phase boundary Eq. (fTTj) and the fact 
that (J z /Ji) x (ojo/t) [(J2/J1) x (u>o/t)] is monotonically 
increasing (decreasing) function of g for g > 1. 

Away from half-filling, the system shows markedly dif- 
ferent behavior compared to the half-filled situation. For 
N p < 4 in the 4x4 lattice considered here, without ac- 
tually presenting the details of the calculations, we first 
note that there is no evidence of a phase transition (for 
both J 2 = and J 2 ^ 0). 

In Fig. g] (c) drawn for N p — 5, although S*(ir, tt) dis- 
plays a CDW transition at a critical value g c — 2.45, n s 
does not go to zero even at large values of g considered. 
Furthermore, we see clearly from Fig. [2] (a) that, above 
this critical value of g, the curvature of the free energy 
curves suggests that the system at N p = 5 is an inhomo- 
geneous mixture of CDW-state and SF-state. Thus away 
from half-filling, at small values of the adiabaticity t/uio, 
our HCB-system undergoes a transition from a SF-state 
to a PS-state at a critical b-p coupling strength [similar 
to the xxz- model in Fig. [T] (b)]. 

However for t/u>o not too small, when NNN hopping is 
present, the system shows a strikingly new behavior for a 
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certain region of the g-parameter space. Let us consider 
the system at N p = 5, t/ui = 1.0, J\ ^ 0, and J2 7^ 0. 
Fig. |4] (d) shows that, above g ss 1.85, the system enters 
a CDW state (as can be seen from the structure factor); 
however, it continues to have a SF character as reflected 
by the finite value of n s . Furthermore, Fig. [3]reveals that 
the curve is concave, i.e., the system is PS, only above 
g = 2.0. This simultaneous presence of CDW and SF 
states, without any inhomogcncity (for 1.85 < g < 2.1), 
implies that the system is a supersolid. Similarly, for 6 
and 7 particles as well, we find that the system undergoes 
transition from a SF-state to a SS-state and then to a PS- 
state. This behavior is displayed in Fig. Q] (a). 

Finally, we shall present the interesting case of Jy = 
and Ji 7^ as a means of understanding the SS phase in 
the phase diagram of Fig. Q] (a). The physical scenario, 
when Ji can be negligibly small compared to J2, has been 
addressed in Ref. [20] for cooperative electron-phonon in- 
teraction in an one-dimensional system. When J\ = 0, 
for large values of nearest-neighbor repulsion, it is quite 
natural that all the particles will occupy a single sub- 
lattice. However, the dramatic jump (at a critical value 
of g), from an equal occupation of both sub-lattices to 
a single sub-lattice occupation, is quite unexpected (see 
Fig. [5|). For a half- filled system, above a critical point, all 
the particles get localized which results in an insulating 
state. This can be seen from Fig. [5] (a). One can see that 
(at g w 1.23) the structure factor dramatically jumps to 
its maximum value, while n s drops to zero and nb takes 
the limiting value of 1/16 for reasons discussed earlier. 
This shows that above g « 1.23, the system is in an in- 
sulating state with one sub-lattice being completely full. 
However, away from half-filling, the system conducts per- 
fectly while occupying a single sub-lattice because of the 
presence of holes in the sub-lattice. For instance, from 
Fig. E] (b) drawn for N p = 5, we see that the structure 
factor jumps to its maximum value at g ~ 1.26, whereas 
n s drops to a finite value which remains constant above 
g = 1.26. We see from Fig. [U(b), based on the curvature 
of the free energy curves, that the 5-particle system does 
not phase separate both above and below the CDW tran- 
sition. In fact, this single-phase-stability is true for any 
filling. This means that the system, at any non-half fill- 
ing and at Ji = 0, exhibits supersolidity above a critical 




value of g\ 

Conclusions: We demonstrated that our BH model 
displays supersolidity in two-dimensions (2D). Our effec- 
tive Hamiltonian of Eq. (JSJ should be realizable in a 
2D optical lattice. Furthermore, our BH model too can 
be mimicked by designing HCB to move as extra parti- 
cles in a 2D array of trapped molecules with HCB cou- 
pled to the breathing mode of the trapped molecules [23| . 
In three-dimensions, supersolidity is more achievable (i) 
in general, due to an increase in the ratio of NNN and 
NN coordination- numbers; and (ii) in particular for bis- 
muthates, because cooperative breathing mode enhances 
the ratio of NNN and NN hoppings. 
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FIG. 5: Comparison of S*(k,tt), nt, and n s when t/u>o = 1.0, 
J 2 / but Ji = and for (a) N p = 8; and (b) N p = 5. 
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